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Summary 

Background The SARS-CoV-2 delta (B.1.617.2) variant is highly transmissible and spreading globally, including in 
populations with high vaccination rates. We aimed to investigate transmission and viral load kinetics in vaccinated 
and unvaccinated individuals with mild delta variant infection in the community. 


Methods Between Sept 13, 2020, and Sept 15, 2021, 602 community contacts (identified via the UK contract-tracing 
system) of 471 UK COVID-19 index cases were recruited to the Assessment of Transmission and Contagiousness of 
COVID-19 in Contacts cohort study and contributed 8145 upper respiratory tract samples from daily sampling for up 
to 20 days. Household and non-household exposed contacts aged 5 years or older were eligible for recruitment if they 
could provide informed consent and agree to self-swabbing of the upper respiratory tract. We analysed transmission 
risk by vaccination status for 231 contacts exposed to 162 epidemiologically linked delta variant-infected index cases. 
We compared viral load trajectories from fully vaccinated individuals with delta infection (n=29) with unvaccinated 
individuals with delta (n=16), alpha (B.1.1.7; n=39), and pre-alpha (n=49) infections. Primary outcomes for the 
epidemiological analysis were to assess the secondary attack rate (SAR) in household contacts stratified by contact 
vaccination status and the index cases' vaccination status. Primary outcomes for the viral load kinetics analysis were 
to detect differences in the peak viral load, viral growth rate, and viral decline rate between participants according to 
SARS-CoV-2 variant and vaccination status. 


Findings The SAR in household contacts exposed to the delta variant was 25% (95% CI 18—33) for fully vaccinated 
individuals compared with 3826 (24—53) in unvaccinated individuals. The median time between second vaccine dose and 
study recruitment in fully vaccinated contacts was longer for infected individuals (median 101 days [IQR 74—120]) than 
for uninfected individuals (64 days [32-97], p=0-001). SAR among household contacts exposed to fully vaccinated index 
cases was similar to household contacts exposed to unvaccinated index cases (25% [95% CI 15—35] for vaccinated vs 23% 
[15—31] for unvaccinated). 12 (39%) of 31 infections in fully vaccinated household contacts arose from fully vaccinated 
epidemiologically linked index cases, further confirmed by genomic and virological analysis in three index case-contact 
pairs. Although peak viral load did not differ by vaccination status or variant type, it increased modestly with age 
(difference of 0-39 [95% credible interval —0. 03 to 0-79] in peak log, viral load per mL between those aged 10 years and 
50 years). Fully vaccinated individuals with delta variant infection had a faster (posterior probability >0-84) mean rate of 
viral load decline (0-95 log, copies per mL per day) than did unvaccinated individuals with pre-alpha (0-69), alpha (0-82), 
or delta (0-79) variant infections. Within individuals, faster viral load growth was correlated with higher peak viral load 
(correlation 0-42 [95% credible interval 0-13 to 0-65]) and slower decline (-0-44 [-0- 67 to —0-18]). 


Interpretation Vaccination reduces the risk of delta variant infection and accelerates viral clearance. Nonetheless, fully 
vaccinated individuals with breakthrough infections have peak viral load similar to unvaccinated cases and can 
efficiently transmit infection in household settings, including to fully vaccinated contacts. Host-virus interactions 
early in infection may shape the entire viral trajectory. 


Funding National Institute for Health Research. 


Copyright © 2021 The Author(s). Published by Elsevier Ltd. This is an Open Access article under the CC BY-NC-ND 
4.0 license. 


Introduction 
While the primary aim of vaccination is to protect 
individuals against severe COVID-19 disease and its 


consequences, the extent to which vaccines reduce 
onward transmission of SARS-CoV-2 is key to containing 
the pandemic. This outcome depends on the ability of 
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Research in context 


Evidence before this study 

The SARS-CoV-2 delta variant is spreading globally, including in 
populations with high vaccination coverage. While vaccination 
remains highly effective at attenuating disease severity and 
preventing death, vaccine effectiveness against infection is 
reduced for delta. Determining the extent of transmission from 
vaccinated delta-infected individuals to their vaccinated 
contacts is a public health priority. Comparing the upper 
respiratory tract (URT) viral load kinetics of delta infections 
with those of other variants gives insight into potential 
mechanisms for its increased transmissibility. We searched 
PubMed and medRxiv for articles published between database 
inception and Sept 20, 2021, using search terms describing 
"SARS-CoV-2, delta variant, viral load, and transmission". 

Two studies longitudinally sampled the URT in vaccinated and 
unvaccinated delta variant-infected individuals to compare viral 
load kinetics. In a retrospective study of a cohort of hospitalised 
patients in Singapore, more rapid viral load decline was found 
in vaccinated individuals than unvaccinated cases. However, the 
unvaccinated cases in this study had moderate-to-severe 
infection, which is known to be associated with prolonged 
shedding. The second study longitudinally sampled 
professional USA sports players. Again, clearance of delta viral 
RNA in vaccinated cases was faster than in unvaccinated cases, 
but only 896 of unvaccinated cases had delta variant infection, 
complicating interpretation. Lastly, a report of a single-source 
nosocomial outbreak of a distinct delta sub-lineage in 
Vietnamese health-care workers plotted viral load kinetics 
(without comparison with unvaccinated delta infections) 

and demonstrated transmission between fully vaccinated 
health-care workers in the nosocomial setting. The findings 
might therefore not be generalisable beyond the particular 
setting and distinct viral sub-lineage investigated. 


Added value of this study 
The majority of SARS-CoV-2 transmission occurs in households, 
buttransmission between fully vaccinated individuals in this 


vaccines to protect against infection and the extent to 
which vaccination reduces the infectiousness of break- 
through infections. 

Vaccination was found to be effective in reducing 
household transmission of the alpha variant (B.1.1.7) by 
40-5096, and infected, vaccinated individuals had 
lower viral load in the upper respiratory tract (URT) 
than infections in unvaccinated individuals? which is 
indicative of reduced infectiousness.** However the 
delta variant (B.1.6172), which is more transmissible than 
the alpha variant,” is now the dominant strain worldwide. 
After a large outbreak in India, the UK was one of the first 
countries to report a sharp rise in delta variant infection. 
Current vaccines remain highly effective at preventing 
admission to hospital and death from delta infection, 
However, vaccine effectiveness against infection is reduced 
for delta, compared with alpha,? and the delta variant 


setting has not been shown to date. To ascertain secondary 
transmission with high sensitivity, we longitudinally followed 
index cases and their contacts (regardless of symptoms) in the 
community early after exposure to the delta variant of 
SARS-CoV-2, performing daily quantitative RT-PCR on URT 
samples for 14-20 days. We found that the secondary attack rate 
in fully vaccinated household contacts was high at 25%, but this 
value was lower than that of unvaccinated contacts (38%). 

Risk of infection increased with time in the 2-3 months since the 
second dose of vaccine. The proportion of infected contacts was 
similar regardless of the index cases' vaccination status. 

We observed transmission of the delta variant between fully 
vaccinated index cases and their fully vaccinated contacts in 
several households, confirmed by whole-genome sequencing. 
Peak viral load did not differ by vaccination status or variant 
type but did increase modestly with age. Vaccinated delta cases 
experienced faster viral load decline than did unvaccinated alpha 
or delta cases. Across study participants, faster viral load growth 
was correlated with higher peak viral load and slower decline, 
suggesting that host-virus interactions early in infection shape 
the entire viral trajectory. Since our findings are derived from 
community household contacts in a real-life setting, they are 
probably generalisable to the general population. 


Implications of all the available evidence 

Although vaccines remain highly effective at preventing severe 
disease and deaths from COVID-19, our findings suggest that 
vaccination is not sufficient to prevent transmission of the 
delta variant in household settings with prolonged exposures. 
Our findings highlight the importance of community studies 
to characterise the epidemiological phenotype of new 
SARS-CoV-2 variants in increasingly highly vaccinated 
populations. Continued public health and social measures 

to curb transmission of the delta variant remain important, 
even in vaccinated individuals. 


continues to cause a high burden of cases even in countries 
with high vaccination coverage. Data are scarce on the risk 
of community transmission of delta from vaccinated 
individuals with mild infections. 

Here, we report data from a UK community-based 
study, the Assessment of Transmission and Conta- 
giousness of COVID-19 in Contacts (ATACCC) study, in 
which ambulatory close contacts of confirmed COVID-19 
cases underwent daily, longitudinal URT sampling, with 
collection of associated clinical and epidemiological 
data. We aimed to quantify household transmission of 
the delta variant and assess the effect of vaccination 
status on contacts risk of infection and index 
cases infectiousness, including (1) households with 
unvaccinated contacts and index cases and (2) house- 
holds with fully vaccinated contacts and fully vaccinated 
index cases. We also compared sequentially sampled 
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URT viral RNA trajectories from individuals with non- 
severe delta, alpha, and pre-alpha SARS-CoV-2 infections 
to infer the effects of SARS-CoV-2 variant status—and, 
for delta infections, vaccination status—on transmission 
potential. 


Methods 

Study design and participants 

ATACCC is an observational longitudinal cohort study of 
community contacts of SARS-CoV-2 cases. Contacts of 
symptomatic PCR-confirmed index cases notified to 
the UK contact-tracing system (National Health Service 
Test and Trace) were asked if they would be willing to 
be contacted by Public Health England to discuss 
participation in the study. All contacts notified within 
5 days of index case symptom onset were selected to be 
contacted within our recruitment capacity. Household 
and non-household contacts aged 5 years or older were 
eligible for recruitment if they could provide written 
informed consent and agree to self-swabbing of the URT. 
Further details on URT sampling are given in the 
appendix (p 13). 

The ATACCC study is separated into two study arms, 
ATACCCI and ATACCC2, which were designed to capture 
different waves of the SARS-CoV-2 pandemic. In 
АТАСССІ, which investigated alpha variant and pre-alpha 
cases in Greater London, only contacts were recruited 
between Sept 13, 2020, and March 13, 2021. ATACCC1 
included a pre-alpha wave (September to November, 2020) 
and an alpha wave (December, 2020, to March, 2021). 
In ATACCC2, the study was relaunched specifically to 
investigate delta variant cases in Greater London and 
Bolton, and both index cases and contacts were recruited 
between May 25, and Sept 15, 2021. Early recruitment was 
focused in West London and Bolton because UK incidence 
of the delta variant was highest in these areas." Based 
on national and regional surveillance data, community 
transmission was moderate-to-high throughout most of 
our recruitment period. 

This study was approved by the Health Research 
Authority. Written informed consent was obtained from all 
participants before enrolment. Parents and caregivers 
gave consent for children. 


Data collection 

Demographic information was collected by the study team 
on enrolment. The date of exposure for non-household 
contacts was obtained from Public Health England. 
COVID-19 vaccination history was determined from the 
UK National Immunisation Management System, general 
practitioner records, and self-reporting by study parti- 
cipants. We defined a participant as unvaccinated if they 
had not received a single dose of a COVID-19 vaccine at 
least 7 days before enrolment, partially vaccinated if they 
had received one vaccine dose at least 7 days before study 
enrolment, and fully vaccinated if they had received 
two doses of a COVID-19 vaccine at least 7 days before 


study enrolment. Previous literature was used to 
determine the 7-day threshold for defining vaccination 
status." We also did sensitivity analyses using a 14-day 
threshold. The time interval between vaccination and 
study recruitment was calculated. We used WHO criteria” 
to define symptomatic status up to the day of study 
recruitment. Symptomatic status for incident cases— 
participants who were PCR-negative at enrolment and 
subsequently tested positive—was defined from the day of 
the first PCR-positive result. 


Laboratory procedures 

SARS-CoV-2 quantitative RT-PCR, conversion of ORFlab 
and envelope (E-gene) cycle threshold values to viral 
genome copies, whole-genome sequencing, and lineage 
assignments are described in the appendix (pp 13-14). 


Outcomes 
Primary outcomes for the epidemiological analysis were 
to assess the secondary attack rate (SAR) in household 
contacts stratified by contact vaccination status and the 
index cases’ vaccination status. Primary outcomes for the 
viral load kinetics analysis were to detect differences in 
the peak viral load, viral growth rate, and viral decline 
rate between participants infected with pre-alpha versus 
alpha versus delta variants and between unvaccinated 
delta-infected participants and vaccinated delta-infected 
participants. 

We assessed vaccine effectiveness and susceptibility to 
SARS-CoV-2 infection stratified by time elapsed since 
receipt of second vaccination as exploratory analyses. 


Statistical analysis 

To model viral kinetics, we used a simple phenomeno- 
logical model of viral titre? during disease pathogenesis. 
Viral kinetic parameters were estimated on a participant- 
specific basis using a Bayesian hierarchical model to fit 
this model to the entire dataset of sequential cycle 
threshold values measured for all participants. For the 
19 participants who were non-household contacts of index 
cases and had a unique date of exposure, the cycle 
threshold data were supplemented by a pseudo-absence 
data point (ie, undetectable virus) on the date of exposure. 
Test accuracy and model misspecification were modelled 
with a mixture model by assuming there was a probability 
p of a test giving an observation drawn from a (normal) 
error distribution and probability 1—р of it being drawn 
from the true distribution. 

The hierarchical structure was represented by grouping 
participants based on the infecting variant and 
their vaccination status. A single-group model was fitted, 
which implicitly assumes that viral kinetic parameters 
vary by individual but not by variant or vaccination 
status. A four-group model was also explored, where 
groups 1, 2, 3, and 4 represent pre-alpha, alpha, 
unvaccinated delta, and fully vaccinated delta, 
respectively. We fitted a correlation matrix between 
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> contacts household contacts 
153 household 
26 non-household 
53 delta-exposed, PCR-positive : 9 with transmission from fully : Total: 
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Figure 1: Recruitment, SARS-CoV-2 infection, variant status, and vaccination history for ATACCC study participants 
(A) Study recruitment and variant status confirmed by whole-genome sequencing (ATACCC1 and ATACCC2 combined). (B) ATACCC2: delta-exposed contacts included in secondary attack rate 
calculation (table 1) and transmission assessment (table 2). NHS=National Health Service. *All index cases were from ATACCC2. ТАП contacts. £The two earliest PCR-positive cases from the 





ATACCC2 cohort (one index case and one contact) were confirmed as having the alpha variant on whole-genome sequencing (recruited on May 28, 2021). This alpha variant-exposed, 

PCR-positive contact is excluded from figure 1B. SOne PCR-negative contact had no vaccination status data available and one PCR-negative contact's index case had no vaccination data available. 
Vaccination data were available for 138 index cases of 163. |[The contacts of these 15 index cases are included within the 232 total contacts. ** These three index cases without contacts are only 
included in the viral load kinetics analysis (figure 3) and are not included in tables 1 and 2. 


participant-specific kinetic parameters to allow us to 
examine whether there is within-group correlation 
between peak viral titre, viral growth rate, and viral 
decline rate. Our initial model selection, using leave-one- 
out cross-validation, selected a four-group hierarchical 
model with fitted correlation coefficients between 
individual-level parameters determining peak viral load 


and viral load growth and decline rates (appendix p 5). 
However, resulting participant-specific estimates of peak 
viral load (but not growth and decline rates) showed a 
marked and significant correlation with age in the 
exploratory analysis, which motivated examination of 
models where mean peak viral load could vary with age. 
The most predictive model overall allowed mean viral 
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load growth and decline rates to vary across the 
four groups, with mean peak viral load common to all 
groups but assumed to vary linearly with the logarithm 
of age (appendix p 5). We present peak viral loads for the 
reference age of 50 years with 95% credible intervals 
(95% CrIs). 50 years was chosen as the reference age as it 
is typical of the ages of the cases in the whole dataset and 
the choice of reference age made no difference in the 
model fits or judgment of differences between the 
groups. 

We computed group-level population means and 
within-sample group means of log peak viral titre, viral 
growth rate, and viral decline rate. Since posterior 
estimates of each of these variables are correlated across 
groups, overlap in the credible intervals of an estimate for 
one group with that for another group does not necessarily 
indicate no significant difference between those groups. 
We, therefore, computed posterior probabilities, pp, 
that these variables were larger for one group than 
another. For our model, Bayes factors can be computed 
as pp/(1-pp). We only report population (group-level) 
posterior probabilities greater than 0-75 (corresponding 
to Bayes factors >3) as indicating at least moderate 
evidence of a difference. 

For vaccine effectiveness, we defined the estimated 
effectiveness at preventing infection, regardless of 
symptoms, with delta in the household setting as 1- SAR 
(fully vaccinated) / SAR (unvaccinated). 


Role of the funding source 

The funder of the study had no role in study design, data 
collection, data analysis, data interpretation, or writing of 
the report. 


Results 

Between Sept 13, 2020, and Sept 15, 2021, 621 community- 
based participants (602 contacts and 19 index cases) from 
471 index notifications were prospectively enrolled in 
the ATACCCI and ATACCC2 studies, and contributed 
8145 URT samples. Of these, ATACCCI enrolled 
369 contacts (arising from 308 index notifications), and 
ATACCC2 enrolled 233 contacts (arising from 163 index 
notifications) and 19 index cases. SARS-CoV-2 RNA was 
detected in 163 (2696) of the 621 participants. Whole- 
genome sequencing of PCR-positive cases confirmed 
that 71 participants had delta variant infection (18 index 
cases and 53 contacts), 42 had alpha variant infection 
(one index case and 41 contacts), and 50 had pre-alpha 
variant infection (all contacts; figure 1A). 

Of 163 PCR-positive participants, 89 (55%) were female 
and 133 (8296) were White. Median age was 36 years 
(IQR 26-50). Sex, age, ethnicity, body-mass index 
(BMI) distribution, and the frequency of comorbidities 
were similar among those with delta, alpha, and 
pre-alpha infection, and for vaccinated and unvaccinated 
delta-infected participants, except for age and sex 
(appendix pp 2-3). There were fewer unvaccinated 
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Total PCR positive PCR negative SAR (95% CI) p value 
Contacts 
All 231 53 178 23 (18-29) NA 
Fully vaccinated 140 31 109 22 (16-30) 0-16 
Unvaccinated 44 15 29 34 (22-49) 
Partially vaccinated 47 7 40 15 (7-28) NA 
Household contacts 
All 205 53 152 26 (20-32) NA 
Fully vaccinated 126 31 95 25 (18-33) 0-17 
Unvaccinated 40 15 25 38 (24-53) 
Partially vaccinated 39 7 32 18 (9-33) NA 
X test was performed to calculate p values for differences in SAR between fully vaccinated and unvaccinated cases. 
One PCR-negative contact who withdrew from the study without vaccination status information was excluded. 
NA-not applicable. SAR=secondary attack rate. 
Table 1: SAR in contacts of delta-exposed index cases recruited to the ATACCC2 study 





females than males (p=0-04) and, as expected from the 
age-prioritisation ofthe UK vaccine roll-out, unvaccinated 
participants infected with the delta variant were 
significantly younger (p<0-001; appendix p 3). Median 
time between exposure to the index case and study 
enrolment was 4 days (IQR 4—5). All participants had 
non-severe ambulatory illness or were asymptomatic. 
The proportion of asymptomatic cases did not differ 
among fully vaccinated, partially vaccinated, and un- 
vaccinated delta groups (appendix p 3). 

No pre-alpha-infected and only опе alpha-infected 
participant had received a COVID-19 vaccine before study 
enrolment. Of 71 delta-infected participants (of whom 
18 were index cases), 23 (3296) were unvaccinated, 
ten (14%) were partially vaccinated, and 38 (54%) were fully 
vaccinated (figure 1A; appendix p 3). Of the 38 fully 
vaccinated delta-infected participants, 14 had received 
the BNT162b2 mRNA vaccine (Pfizer-BioNTech), 23 the 
ChAdOx1 nCoV-19 adenovirus vector vaccine (Oxford- 
AstraZeneca), and one the CoronaVac inactivated whole- 
virion vaccine (Sinovac). 

Itis highly probable thatall but one ofthe 233 ATACCC2 
contacts were exposed to the delta variant because they 
were recruited when the regional prevalence of delta was 
at least 90%, and mostly 95-99% (figure 1B). Of these, 
206 (89%) were household contacts (in 127 households), 
and 26 (11%) were non-household contacts. Distributions 
of age, ethnicity, BMI, smoking status, and comorbidities 
were similar between PCR-positive and PCR-negative 
contacts (appendix p 4). The median time between 
second vaccine dose and study recruitment in fully 
vaccinated contacts with delta variant infection was 
74 days (IQR 35-105; range 16-201), and this was 
significantly longer in PCR-positive contacts than in 
PCR-negative contacts (101 days [IQR 74—120] vs 64 days 
[32-97], respectively, p=0-001; appendix p 4). All 
53 PCR-positive contacts were exposed in household 
settings and the SAR for all delta variant-exposed 
household contacts was 26% (95% CI 20-32). SAR was 
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contacts exposed to each category of index case. 


contacts (n=204)* (п=125) (n=39) (n=40) 
PCR positive PCR negative PCR positive PCR negative PCR positive PCR negative 
(n=31) (n=94) (n=7) (n=32) (n=15) (n=25) 
Fully vaccinated index cases (n=50) 69 12) 31 1 8 4 13 
Partially vaccinated index cases (п=25) 35 7 12 3 10 3 0 
Unvaccinated index cases (n=63) 100 12 51 3 14 8 12 


Non-household exposed contacts (n=24, all PCR negative) were excluded. One PCR-negative household contact who withdrew from the study without vaccination status 
information was excluded. One PCR-negative household contact who could not be linked to their index case was also excluded. *The rows below show the number of 








Table 2: Comparison of vaccination status of the 138 epidemiologically linked PCR-positive index cases for 204 delta variant-exposed household contacts 








not significantly higher in unvaccinated (38%, 95% CI 
24-53) than fully vaccinated (25%, 18-33) household 
contacts (table 1). We estimated vaccine effectiveness at 
preventing infection (regardless of symptoms) with delta 
in the household setting to be 34% (bootstrap 95% CI 
—15 to 60). Sensitivity analyses using a 14 day threshold 
for time since second vaccination to study recruitment to 
denote fully vaccinated did not materially affect our 
estimates of vaccine effectiveness or SAR (data not 
shown). Although precision is restricted by the small 
sample size, this estimate is broadly consistent with 
vaccine effectiveness estimates for delta variant infection 
based on larger datasets.?'97 

The vaccination status of 138 epidemiologically linked 
index cases of 204 delta variant-exposed household 
contacts was available (figure 1B, table 2). The SAR in 
household contacts exposed to fully vaccinated index 
cases was 25% (95% CI 15—35; 17 of 69), which is similar 
to the SAR in household contacts exposed to unvaccinated 
index cases (23% [15-31]; 23 of 100; table 2). The 
53 PCR-positive contacts arose from household exposure 
to 39 PCR-positive index cases. Of these index cases who 
gave rise to secondary transmission, the proportion who 
were fully vaccinated (15 [38%] of 39) was similar to the 
proportion who were unvaccinated (16 [41%] of 39). The 
median number of days from the index cases’ second 
vaccination to the day of recruitment for their respective 
contacts was 73 days (IQR 38-116). Time interval did not 
differ between index cases who transmitted infection to 
their contacts and those who did not (94 days [IQR 62-112] 
and 63 days [35-117], respectively; p=0-43). 

18 of the 163 delta variant-infected index cases that led 
to contact enrolment were themselves recruited to 
ATACCC2 and serial URT samples were collected from 
them, allowing for more detailed virology and genome 
analyses. For 15 ofthese, their contacts were also recruited 
(13 household contacts and two non-household contacts). 
A corresponding PCR-positive household contact was 
identified for four of these 15 index cases (figure 1B). 
Genomic analysis showed that index-contact pairs were 
infected with the same delta variant sub-lineage in 
these instances, with one exception (figure 2A). In 
one household (number 4), an unvaccinated index case 
transmitted the delta variant to an unvaccinated contact, 


while another partially vaccinated contact was infected 
with a different delta sub-lineage (which was probably 
acquired outside the household). In the other three 
households (numbers 1-3), fully vaccinated index cases 
transmitted the delta variant to fully vaccinated household 
contacts, with high viral load in all cases, and temporal 
relationships between the viral load kinetics that were 
consistent with transmission from the index cases to 
their respective contacts (figure 2B). 

Inclusion criteria for the modelling analysis selected 
133 participant's viral load RNA trajectories from 
163 PCR-positive participants (49 with the pre-alpha 
variant, 39 alpha, and 45 delta; appendix p 14). Of the 
45 delta cases, 29 were fully vaccinated and 16 were 
unvaccinated; partially vaccinated cases were excluded. 
Of the 133 included cases, 29 (22%) were incident 
(ie, PCR negative at enrolment converting to PCR positive 
subsequently) and 104 (78%) were prevalent (ie, already 
PCR positive at enrolment). 15 of the prevalent cases had 
a clearly resolvable peak viral load. Figure 3 shows 
modelled viral RNA (ORF1ab) trajectories together with 
the viral RNA copy numbers measured for individual 
participants. The E-gene equivalent is shown in the 
appendix (p 2). Estimates derived from E-gene cycle 
threshold value data (appendix pp 5, 7, 9, 11) were similar 
to those for ORF lab. 

Although viral kinetics appear visually similar for all four 
groups of cases, we found quantitative differences in 
estimated viral growth rates and decline rates (tables 3, 4). 
Population (group-level) estimates of mean viral load 
decline rates based on ORF1ab cycle threshold value data 
varied in the range of 0-69-0-95 log, units per mL 
per daxes 4; appendix p 10), indicating that a typical 
10-day period was required for viral load to decline from 
peak to undetectable. A faster decline was seen in the alpha 
(pp=0-93), unvaccinated delta (pp=0-79), and fully 
vaccinated delta (pp=0-99) groups than in the pre-alpha 
group. The mean viral load decline rate of the fully 
vaccinated delta group was also faster than those of the 
alpha group (pp=0-84) and the unvaccinated delta group 
(pp=0-85). The differences in decline rates translate into a 
difference of about 3 days in the mean duration of the 
decline phase between the pre-alpha and delta vaccinated 


groups. 
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Figure 2: Virological, epidemiological, and genomic evidence for transmission of the SARS-CoV-2 delta variant (B.1.617.2) in households 

(A) Genomic analysis of the four households with lineage-defining mutations for delta? and additional mutations within ORFs displayed to give insight into whether 
strains from individuals within the household are closely related. Lineages AY.4 and AY.9 are sub-lineages of delta. (B) Viral trajectories and vaccination status of the 
four index cases infected with the delta variant for whom infection was detected in their epidemiologically linked household contacts. All individuals had non-severe 
disease. Each plot shows an index case and their household contacts. Undetectable viral load measurements are plotted at the limit of detection (10*^). C=contact. 


I-index case. FV=fully vaccinated. ORF=open reading frame. PV=partially vaccinated. 


Viral load growth rates were substantially faster than 
decline rates, varying in the range of 2-69-3-24 log, 
units per mL per day between groups, indicating that a 
typical 3-day period was required for viral load to 


U-unvaccinated. 


grow from undetectable to peak. Our power to infer 
differences in growth rates between groups was more 
restricted than for viral decline, but there was moderate 
evidence (pp=0-79) that growth rates were lower for 
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VL growth rate 


(95% Crl), од 
units per day 


Posterior probability 
estimate is less than 
pre-alpha 


Posterior probability 
estimate is less than 
alpha 


Posterior probability 
estimate is less than 
delta (unvaccinated) 


Posterior probability 
estimate is less than 
delta (fully 


vaccinated) 





Pre-alpha (n=49) 
Alpha (n=39) 


Delta, unvaccinated (n=16) 


3:24 (1:78-6-14) 
3-13 (1:76-5:94) 0:56 
2:81 (1:47-5:47) 073 


Delta, fully vaccinated (n=29) 2-69 (1:51-5:17) 0-79 


errors of <0-01. VL=viral load. Crl=credible interval. 


VL growth rates are shown as within-sample posterior mean estimates. Remaining columns show population (group-level) posterior probabilities that the estimate on that 
row is less than an estimate for a different group. Posterior probabilities are derived from 20 000 posterior samples and have sampling 


0-44 0-27 0-21 

n 0-32 0-25 
0-68 E 0-44 
0-75 0-56 





Table 3: Estimates of VL growth rates for pre-alpha, alpha, and delta (unvaccinated and fully vaccinated) cases, derived from ORF1ab cycle threshold data 








VL decline rate 
(95% Сп), од 


Posterior probability Posterior probability Posterior probability Posterior probability 
estimate is larger 


estimate is larger estimate is larger estimate is larger 





errors of <0-01. VL=viral load. Crl=credible interval. 


units per day than pre-alpha than alpha than delta than delta (fully 
(unvaccinated) vaccinated) 
Pre-alpha (n=49) 0-69 (0-58-0-81) 0-07 0-21 0-01 
Alpha (n=39) 0-82 (0-67-1-01) 0-93 0-60 0-16 
Delta, unvaccinated (n=16) 0-79 (0-59-1-04) 0-79 0-40 - 0-15 
Delta, fully vaccinated (n=29) 0-95 (0:76-1:18) 0-99 0-84 0-85 


VL decline rates are shown as within-sample posterior mean estimates. Remaining columns show population (group-level) posterior probabilities that the estimate on that 
row is less than an estimate for a different group. Posterior probabilities are derived from 20 000 posterior samples and have sampling 








Table 4: Estimates of VL decline rates for pre-alpha, alpha, and delta (unvaccinated and fully vaccinated) cases, derived from ORF1ab cycle threshold data 








those in the vaccinated delta group than in the pre-alpha 
group. 

We estimated mean peak viral load for 50-year-old 
adults to be 8-14 (95% С 7-95 to 8-32) log, copies 
per mL, but peak viral load did not differ by variant or 
vaccination status. However, we estimated that peak viral 
load increases with age (pp=0-96 that the slope of peak 
viral load with log[age] was >0), with an estimated 
slope of 0-24 (95% Cr1 -0-02 to 0-49) log, copies per mL 
per unit change in log(age). This estimate translates to a 
difference of 0-39 (-0-03 to 0-79) in mean peak log, 
copies per mL between those aged 10 years and 50 years. 

Within-group individual participant estimates of viral 
load growth rate were positively correlated with peak viral 
load, with a correlation coefficient estimate of 0-42 
(95% CrI 0-13 to 0-65; appendix p 8). Hence, individuals 
with faster viral load growth tend to have higher peak 
viral load. The decline rate of viral load was also negatively 
correlated with viral load growth rate, with a correlation 
coefficient estimate of —0-44 (95% CrI —0-67 to —0-18), 
illustrating that individuals with faster viral load growth 
tend to experience slower viral load decline. 


Discussion 

Households are the site of most SARS-CoV-2 transmission 
globally.’ In our cohort of densely sampled household 
contacts exposed to the delta variant, SAR was 3896 in 
unvaccinated contacts and 25% in fully vaccinated 
contacts. This finding is consistent with the known 
protective effect of COVID-19 vaccination against 


infection.^ Notwithstanding, these findings indicate 
continued risk of infection in household contacts despite 
vaccination. Our estimate of SAR is higher than that 
reported in fully vaccinated household contacts exposed 
before the emergence of the delta variant."" The time 
interval between vaccination and study recruitment was 
significantly higher in fully vaccinated PCR-positive 
contacts than fully vaccinated PCR-negative contacts, 
suggesting that susceptibility to infection increases with 
time as soon as 2-3 months after vaccination—consistent 
with waning protective immunity. This potentially 
important observation is consistent with recent large-scale 
data and requires further investigation." Household SAR 
for delta infection, regardless of vaccination status, 
was 2696 (9596 CI 20-32), which is higher than estimates 
of UK national surveillance data (10.896 [10-7-10-9]).” 
However, we sampled contacts daily regardless of 
symptomatology, to actively identify infection with 
high sensitivity. By contrast, symptom-based, single- 
timepoint surveillance testing probably underestimates 
the true SAR, and potentially also overestimates vaccine 
effectiveness against infection. 

We identified similar SAR (25%) in household contacts 
exposed to fully vaccinated index cases as in those exposed 
to unvaccinated index cases (23%). This finding indicates 
that breakthrough infections in fully vaccinated people can 
efficiently transmit infection in the household setting. We 
identified 12 household transmission events between fully 
vaccinated index case-contact pairs; for three of these, 
genomic sequencing confirmed that the index case and 
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contact were infected by the same delta variant sub-lineage, 
thus substantiating epidemiological data and temporal 
relationships of viral load kinetics to provide definitive 
evidence for secondary transmission. To our knowledge, 
one other study has reported that transmission of the delta 
variant between fully vaccinated people was a point-source 
nosocomial outbreak—a single health-care worker with a 
particular delta variant sub-lineage in Vietnam.” 

Daily longitudinal sampling of cases from early (median 
4 days) after exposure for up to 20 days allowed us to 
generate high-resolution trajectories of URT viral load over 
the course of infection. To date, two studies have sequen- 
tially sampled community cases of mild SARS-CoV-2 
infection, and these were from highly specific population 
groups identified through asymptomatic screening 
programmes (eg, for university staff and students? and 
for professional athletes”). 

Our most predictive model of viral load kinetics 
estimated mean peak log,, viral load per mL of 8-14 
(95% CrI 7-95-8-32) for adults aged 50 years, which is 
very similar to the estimate from a 2021 study using 
routine surveillance data.” We found no evidence of 
variation in peak viral load by variant or vaccination 
status, but we report some evidence of modest but 
significant (pp=0-95) increases in peak viral load with 
age. Previous studies of viral load in children and 
adults*^"* have not used such dense sequential sampling 
of viral load and have, therefore, been restricted in their 
power to resolve age-related differences; the largest such 
study” reported a similar difference between children 
and adults to the one we estimated. We found the rate of 
viral load decline was faster for vaccinated individuals 
with delta infection than all other groups, and was faster 
for individuals in the alpha and unvaccinated delta 
groups than those with pre-alpha infection. 

For all variant vaccination groups, the variation 
between participants seen in viral load kinetic parameter 
estimates was substantially larger than the variation in 
mean parameters estimated between groups. The 
modest scale of differences in viral kinetics between 
fully vaccinated and unvaccinated individuals with 
delta infection might explain the relatively high rates of 
transmission seen from vaccinated delta index cases in 
our study. We found no evidence of lower SARs from 
fully vaccinated delta index cases than from unvaccinated 
ones. However, given that index cases were identified 
through routine symptomatic surveillance, there might 
have been a selection bias towards identifying untypically 
symptomatic vaccine breakthrough index cases. 

The differences in viral kinetics we found between the 
pre-alpha, alpha, and delta variant groups suggest some 
incremental, but potentially adaptive, changes in viral 
dynamics associated with the evolution of SARS-CoV-2 
towards more rapid viral clearance. Our study provides 
the first evidence that, within each variant or vaccination 
group, viral growth rate is positively correlated with peak 
viral load, but is negatively correlated with viral decline 


rate. This finding suggests that individual infections 
during which viral replication is initially fastest generate 
the highest peak viral load and see the slowest viral 
clearance, with the latter not just being due to the higher 
peak. Mechanistically, these data suggest that the host and 
viral factors determining the initial growth rate of 
SARS-CoV-2 have a fundamental effect on the trajectory 
throughout infection, with faster replication being more 
difficult (in terms of both peak viral load and the 
subsequent decline of viral load) for the immune response 
to control. Analysis of sequentially sampled immune 
markers during infection might give insight into the 
immune correlates of these early differences in infection 
kinetics. It is also possible that individuals with the 
fastest viral load growth and highest peaks contribute 
disproportionately to community transmission, a hypo- 
thesis that should be tested in future studies. 

Several population-level, single-timepoint sampling 
studies using routinely available data have found no major 
differences in cycle threshold values between vaccinated 
and unvaccinated individuals with delta variantinfection.°”* 
However, as the timepoint of sampling in the viral trajectory 
is unknown, this restricts the interpretation of such results. 
Two other studies longitudinally sampled vaccinated and 
unvaccinated individuals with delta variant infection.” 
A retrospective cohort of hospitalised patients in Singapore” 
also described a faster rate of viral decline in vaccinated 
versus unvaccinated individuals with delta variant, reporting 
somewhat larger differences in decline rates than we 
estimated here. However, this disparity might be accounted 
for by the higher severity of illness in unvaccinated 
individuals in the Singaporean study (almost two-thirds 
having pneumonia, one-third requiring COVID-19 treat- 
ment, and a fifth needing oxygen) than in our study, given 
that longer viral shedding has been reported in patients 
with more severe illness.” A longitudinal sampling 
study in the USA reported that pre-alpha, alpha, and 
delta variant infections had similar viral trajectories.^ The 
study also compared trajectories in vaccinated and 
unvaccinated individuals, reporting similar proliferation 
phases and peak cycle threshold values, but more rapid 
clearance of virus in vaccinated individuals. However, this 
study in the USA stratified by vaccination status and variant 
separately, rather than jointly, meaning vaccinated 
individuals with delta infection were being compared with, 
predominantly, unvaccinated individuals with pre-alpha 
and alpha infection. Moreover, sampling was done as part of 
a professional sports player occupational health screening 
programme, making the results not necessarily repre- 
sentative of typical community infections. 

Our study has limitations. First, we recruited only 
contacts of symptomatic index cases as our study 
recruitment is derived from routine contact-tracing 
notifications. Second, index cases were defined as the first 
household member to have a PCR-positive swab, but we 
cannot exclude the possibility that another household 
member might already have been infected and transmitted 
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to the index case. Third, recording of viral load trajectories 
is subject to left censoring, where the growth phase in 
prevalent contacts (already PCR-positive at enrolment) was 
missed for a proportion of participants. However, we 
captured 29 incident cases and 15 additional cases on the 
upslope of the viral trajectory, providing valuable, 
informative data on viral growth rates and peak viral load 
in a subset of participants. Fourth, owing to the 
age-stratified rollout of the UK vaccination programme, 
the age of the unvaccinated, delta variant-infected parti- 
cipants was lower than that of vaccinated participants. 
Thus, age might be a confounding factor in our results 
and, as discussed, peak viral load was associated with age. 
However, it is unlikely that the higher SAR observed in the 
unvaccinated contacts would have been driven by younger 
age rather than the absence of vaccination and, to our 
knowledge, there is no published evidence showing 
increased susceptibility to SARS-CoV-2 infection with 
decreasing age.” Finally, although we did not perform viral 
culture here—which is a better proxy for infectiousness 
than RT-PCR—two other studies”? have shown cultivable 
virus from around two-thirds of vaccinated individuals 
infected with the delta variant, consistent with our 
conclusions that vaccinated individuals still have the 
potential to infect others, particularly early after infection 
when viral loads are high and most transmission is 
thought to occur.” 

Our findings help to explain how and why the 
delta variant is being transmitted so effectively in 
populations with high vaccine coverage. Although 
current vaccines remain effective at preventing severe 
disease and deaths from COVID-19, our findings suggest 
that vaccination alone is not sufficient to prevent all 
transmission of the delta variant in the household 
setting, where exposure is close and prolonged. 
Increasing population immunity via booster programmes 
and vaccination of teenagers will help to increase the 
currently limited effect of vaccination on transmission, 
but our analysis suggests that direct protection of 
individuals at risk of severe outcomes, via vaccination 
and non-pharmacological interventions, will remain 
central to containing the burden of disease caused by the 
delta variant. 
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